Isospin Dependence of the Spin Orbit Force and Effective Nuclear Potentials 

M.M. SHARMA\ G. LALAZISSIS^, J. KONIG^, and P. 

RING^ 

^Max Planck Institut fiir Astrophysik, Karl-Schwarzschildstrasse 1, D-85740 Garching, Germany 
Physik- Department der Technischen Universitdt Miinchen, D-85747 Garching, Germany 

(Received 16 October, 1993) 

The isospin dependence of the spin-orbit potential is investigated for an effective 
Skyrme-hke energy functional suitable for density dependent Hartree-Fock calcula- 
tions. The magnitude of the isopin dependence is obtained from a fit to experimental 
data on finite spherical nuclei. It is found to be close to that of relativistic Hartree 
models. Consequently, the anomalous kink in the isotope shifts of Pb nuclei is well 
reproduced. 
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The Hartree-Fock approach based upon phenomenological density dependent forces 
has proved to be very successful in the microscopic description of ground state properties 
of nuclear matter and of finite nuclei over the entire periodic table. In all these calculations 
the spin-orbit potential has been assumed to be isospin independent. Its only parameter, 
the strength, is usually adjusted to the experimental spin-orbit splitting in spherical nuclei 
like ^^O. The exchange term, however, causes for nuclei with neutron excess a strong isospin 
dependence of the correponding single-particle spin-orbit field. 

In recent years Relativistic Mean Field (RMF) theory [^] with nonlinear self-interactions 
between the mesons has gained considerable interest for the investigations of low-energy 
phenomena in nuclear structure. With only a few phenomenological parameters such the- 
ories are able to give a quantitative description of ground state properties of spherical and 
deformed nuclei P,^] at and away from the stability line. In addition, excellent agreement 
with experimental data has been found recently also for collective excitations such as giant 
resonances [0] and for twin bands in rotating superdeformed nuclei p] . In many respects the 
relativistic mean-field theory is regarded as similar to the density-dependent Hartree-Fock 
theory of the Skyrme type P,pi| 

Recently, however, detailed investigations of high precision data on nuclear charge radii in 



Pb isotopes [IT^,!!!] and of shell effects at the neutron drip line |113||l^ have shown considerable 
differences between the Skyrme approach and the relativistic mean field theory. In fact, 
density dependent Hartree-Fock calculations with Skyrme [|TT| or Gogny forces [|15[, which 
used so far antisymmetrized, isospin independent spin-orbit interactions, were not able to 
reproduce the kink in the isotope shifts of Pb nuclei (see Fig. p. On the other hand. 



this kink is obtained in the RMF theory without any new adjustment of parameters [jT2 
Another considerable difference has been found in theoretical investigation of shell effects in 
very exotic Zr-isotopes near the neutron drip line: in conventional non-relativistic Skyrme 
calculations the shell gap at isotope ^^^Zr with the magic neutron configuration N = 82 is 
totally smeared out |]T3, whereas relativistic calculations using various parameter sets show 



at N=82 a clear kink in the binding energy as a function of the neutron number [Q. This 



difference is caused by the different spin-orbit splitting of the single particle levels in these 
nuclei. Mass calculations within the Finite Range Droplet Model (FRDM) ||T6[, which are 
based on an isospin dependent spin-orbit term carefully adjusted to experimental data, are 
in excellent agreement with the relativistic predictions. 

This gives us a hint, that one does not really need full relativistic calculations in order to 
understand these differences between conventional density dependent Hartree-Fock calcula- 
tions and RMF theory. A spin-orbit term with a properly chosen isospin dependence might 
represent the essential part of a relativistic calculation. In this letter, we therefore explore 
the isospin dependence of the spin-orbit term in non-relativistic Skyrme calculations and 
analyze the consequences of this on the nuclear properties. We start from the Skyrme type 
force 

\/(l,2)=to(l+a;oP")<5(ri-r2) 

+ h{l + XiP''){6{ri-T2)k^ + h.c.) 

+ t2(l + a;2P")k5(ri -r2)k (1) 

+ it3(l + a;3P'^)5(ri-r2)p" 
o 

+ Wo{l + x^P")(o-« + o-(2))k X 5(ri - r2)k 

with k = |(pi — P2)- In contrast to the conventional Skyrme ansatz, where the energy 
functional contains a Hartree- and a Fock-contribution, we here neglect the exchange (Fock- 
) term for the spin-orbit potential in the last line of Eq. (|l]). Otherwise, the operator P"^ 
would be equivalent to +1 for spin saturated systems. For the rest of the potential the Fock 
terms are included. The associated 11 parameters U, Xi {i = . . .3), Wq, and a of the 
Modified Skyrme Ansatz (MSkA) in Eq. (|I]) are determined by a fit to experimental data 
of finite spherical nuclei. The nuclear properties taken into consideration are the empirical 
binding energies and charge radii of the closed-shell nuclei ^^O, ^°Ca, ^''Zr, and ^°*Pb. In 
order to take into account the variation in isospin we have also included Sn-isotopes ^^^Sn, 
^^''Sn, and the doubly closed nucleus ^'^^Sn as well as one of the lead isotopes ^^^Pb. The 
resulting force and its parameters are presented in Table |. 
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The spin-orbit term in the single-particle field derived from this force has the form 
W^(r)(p X cr) with 

W,(r) = W^Vpr + ,W2Vpr'^r, (2) 

where Pr is the density for neutrons or protons (r = n or p) and Wi = Wo{l + Xw)/2, 
W2 = Wo/2. Conventional Skyrme calculations use a spin-orbit potential without isospin 
dependence and include the Fock term. This leads to x^, = 1 and the relationship W1/W2 = 
2, which is nearly by a factor 2 different from the value 1.0005 obtained within the modified 
Skyrme Ansatz MSkA (see Table |). It is interesting to note that the fit leads to a value of 
Xw very close to zero, which corresponds to the one without the Fock term. 

In order to study the spin-orbit term in the RMF theory we start from the standard 
Lagrangian density 

£ = 7/> {'y{id - g^uj - g^pf - eA) - m - g^a) ^ 
+ l{daf -U{a) --n.^n^'^ + \mlu^' 

+ \my - (3) 

which contains nucleons with mass m. a-, tu-, p-mesons, the electromagnetic field and 
nonlinear self-interactions U{a) of the a-field, 

U{a) = \mla' + \g2a' + \g,a\ (4) 



In a non-relativistic approximation of the corresponding Dirac equation [|T^ we obtain 
the following single particle spin-orbit term: 

^r{r) = ^^{Cl + Cl + Cl)Vp^ 

+ ^^(C' + Cl- C^^^Pr'+r. (5) 

with Cf = {mgi/niiY for i = a,uj,p, which is similar in form to the spin-orbit field (0), but 
which contains r-dependent parameters 
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^^ = ^^^^^^ + ^^ + (6) 
W2 = ^^,{Cl + Cl-Cl). (7) 

The r- dependence drops out in the ratio W1/W2 which is 1.13 for the parameter set NLl 
and 1.10 for the parameter set NL-SH This is only shghtly higher than the value 1.0005 
obtained in the Modified Skyrme Ansatz by the fit to the empirical data (see Table |). The 
absolute size of the spin-orbit term turns out to be r-dependent, which stems from the r- 
dependence of the effective mass m*(r). It can be approximated by m/m* ^ 1 + p{r). 
A more careful consideration would therefore require an explicitly density dependent spin- 
orbit term. It has not been included in the present investigation. 

In Table |I| we show nuclear matter results obtained in the MSkA and compare it with 
the values from the conventional Skyrme force SkM*. The saturation density is obtained 
as po = 0.1531 fm~^ in MSkA. This is the same as the value obtained from an extensive 
fit of the mass formula FRDM [|T^. The binding energy per particle E/A is 16.006 MeV, 
which is close to that of other Skyrme forces. The compression modulus K = 319 MeV 
is somewhat higher than that of the presently adopted Skyrme forces, but it is in good 



agreement with a recent analysis based upon the breathing mode energies ||T^. The third 



derivative e'" = d^{E/A)/dp^ has a value, which is close to that obtained in Ref. |18]. The 



asymmetry energy J is close to the empirical value of 33 MeV. The effective mass m* is in 
good agreement with that of SkM*. 

In Table |ITl| we show binding energies and charge radii obtained in the MSkA for a 
number of spherical nuclei and compare them with those from SkM*. A comparison with 
the empirical values shows, that the binding energies obtained with MSkA have improved 
over those of SkM*. The slightly reduced charge radii in MSkA seem to be connected with 
a slightly higher binding energy, which improves the results for the lighter Pb-isotopes. In 
general the charge radii are improved including that of ^^O. 

In Fig. 1 we show the isotope shifts of Pb nuclei for the Modified Skyrme Ansatz MSkA 
together with experimental data and results of conventional Skyrme calculations. For MSkA 



we observe a clear kink at the double magic nucleus ^°*Pb, whereas the conventional Skyrme 
force SkM* with an isospin independent spin-orbit term gives an almost straight line. For 
the lighter isotopes both theories give excellent agreement, on the heavier side MSkA comes 
closer to the experimental isotope shifts. It may be recalled, that the RMF theory |]12[ is 
successful in reproducing the full size of this kink. MSkA uses a density independent spin- 
orbit force. In contrast the spin-orbit term derived from the RMF theory (see Eq. ^ is 
implicitly density dependent through the density dependence of the effective mass m*(r). A 
density dependence of the spin-orbit term in Skyrme theory might improve the charge radii 
of the heavier Pb-isotopes. 

In Fig. 2 we present binding energies of Zr-nuclei about the neutron drip line as a 
function of the mass number. It is observed that in agreement with earlier investigations 
||r3| within the conventional Skyrme theory shell effects about the closed shell nucleus ^^^Zr 
are weakened considerably. In contrast, RMF-calculations exhibit strong shell effects in the 
drip line region. It has been surmised in Ref. ||I4[ that this is caused by the differences in 
the spin-orbit terms in the two approaches. It is gratifying to see that introduction of an 
isospin dependent (not antisymmetrized) spin-orbit term in Skyrme theory leads to stronger 
shell effects. This is also in agreement with the predictions of the FRDM [|r^, where also 



an isospin dependent spin-orbit term is used. A more quantitative analysis shows that the 
ratio WxlW2 in FRDM is close to 1.06 for the Pb-nuclei and 1.09 for ^'^'^Zx. 

Summarizing, we conclude that the isospin dependence of the spin-orbit term has an 
essential influence on the details of anomalous isotope shifts of Pb-nuclei. A new Modified 
Skyrme Ansatz has been proposed, and the isospin dependence of the spin-orbit strength 
has been determined. The magnitude of this isospin dependence is in agreement with the 
deductions from the relativistic mean-field theory. A reasonably good agreement with the 
experimental data on the binding energies and charge radii has been obtained. The kink in 
the isotope shifts of Pb-nuclei has been obtained in the modified Skyrme ansatz. However, 
the agreement with the empirical isotope shifts for heavy Pb-nuclei is not so good as in 
the RMF theory. This calls for further investigations including a density dependence of the 
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spin-orbit term. 

Finally a remark about the isospin dependence of the spin-orbit term: our investigations 
lead to the interesting result, that the parameter in the Eq. (1) is close to zero. This leads 
practically to an isospin independent single-particle spin-orbit field in Eq. (2). This is in 
agreement with relativistic calculations, where the entire isospin-dependence of the spin-orbit 
field is caused by the parameter Cp in Eq. (5), which is in fact rather small. In contrast the 
two-body spin-orbit potiential in conventional Skyrme theory is isospin independent. The 
exchange term, however, causes a strong isospin dependence in the corresponding single- 
particle spin-orbit field. On the other hand, in the RMF theory the spin-orbit field has its 
origin in Lorentz covariance. There is no contribution from a two-body spin-orbit potential 
and an exchange term is therefore excluded. 

One of the authors (G.A.L) acknowledges support from the E.U., HCM program, con- 
tract: EG/ERB CHBICT-930651 , This work is also supported in part by the Bundesmin- 
isterium fiir Forschung und Technologic under the project 06 TM 733. 
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FIGURES 

FIG. 1. Isotopic shifts in the charge radii of Pb isotopes normalized to the nucleus ^°^Pb as a 
function of the mass number A in the Modified Skyrme Ansatz (MSkA) as compared to the Skyrme 
force SkM* and the empirical data. 



FIG. 2. Binding energy of Zr isotopes at neutron drip line in the MSkA and SkM* calculation. 
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TABLES 



TABLE I. Parameters of the interaction in the modified Skyrme Ansatz (MSkA). 

to = -1200.074 (MeV fm^) xq = 0.187 

h = 396.302 (MeV fm^) xi = 0.018 

t2 = -105.579 (MeV fm^) X2 = -0.059 

t3 = 10631.527 (MeV im^+^°') xs = 0.046 

Wo = 316.38 (MeV fm^) x^ = 0.0005 
a = 0.7557 

TABLE II. Nuclear matter properties obtained in the modified Skyrme Ansatz (MSkA). 

MSkA SkM* 

po 0.1531 fm-3 0.1603 fm'^ 

(E/A)^ 16.006 MeV 15.776 MeV 

K 319.4 MeV 216.7 MeV 

pie'" 100.2 MeV 57.9 MeV 

J 30.0 MeV 30.0 MeV 

m*/m 0.76 0.79 
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TABLE III. The binding energies and charge radii obtained with the Modified Skyrme Ansatz 
(MSkA) in the Hartree-Fock approximation as compared with the normal Skyrme force SkM* and 
the empirical values. 



Nuclei 


Binding Energies 


(MeV) 


Charj 


^e Radii (fm) 




expt. 


MSkA 


SkM* 


expt. 


MSkA 


SkM* 


160 


-127.6 


-128.1 


-127.7 


2.730 


2.742 


2.811 


40Ca 


-342.1 


-342.8 


-341.1 


3.450 


3.468 


3.518 


^«Ca 


-416.0 


-416.8 


-420.1 


3.500 


3.506 


3.537 


90Zr 


-783.9 


-781.9 


-783.0 


4.270 


4.274 


4.296 




-988.7 


-984.0 


-983.4 


4.626 


4.623 


4.619 


i24Sn 


-1050.0 


-1047.9 


-1049.0 


4.673 


4.677 


4.678 


132 Sn 


-1102.9 


-1106.3 


-1110.7 


- 


4.728 


4.727 


200 p]^ 


-1576.4 


-1570.9 


-1568.4 


5.464 


5.465 


5.468 


202pb 


-1592.2 


-1588.3 


-1586.0 


5.473 


5.475 


5.478 


204pb 


-1607.5 


-1605.4 


-1603.4 


5.483 


5.486 


5.489 


206pb 


-1622.3 


-1622.3 


-1620.3 


5.492 


5.496 


5.501 


208pb 


-1636.5 


-1637.8 


-1636.4 


5.503 


5.506 


5.510 


210pb 


-1645.6 


-1645.8 


-1645.6 


5.522 


5.522 


5.520 


212pb 


-1654.5 


-1653.8 


-1654.5 


5.540 


5.537 


5.531 


214pb 


-1663.3 


-1661.6 


-1663.1 


5.558 


5.552 


5.541 
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